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Abstract 


The  design  and  testing  of  an  all-thermistor  (no  thermopile)  pyrgeometer  (LWT)  is  described. 
After  calibration  by  comparison  with  a  Kipp  &  Zonen  CG4,  9.2  months  of  data  show  mean 
differences  of  order  1-2  W  m'2  with  standard  deviations  of  order  7-8  W  m'2.  Approximately  half 
of  the  mean  difference  and  the  standard  deviation  derived  from  40  occasions  when  the  LWT 
readings  were  anomalously  high  for  periods  of  2-14  hours,  principally  at  night.  No  reason  has 
been  found  for  the  anomalous  behavior.  During  the  9.2  months  of  data,  there  were  also  11 
periods  of  a  few  hours  each  when  the  Eppley  PER  indicated  noticeably  higher  flux  values  than 
did  the  Kipp  &  Zonen  CG4. 

The  conclusion  reached  is  that  contemporary  thermistors  allow  temperature  measurements  of 
sufficient  accuracy,  and  the  thermopile  can  be  eliminated  from  pyrgeometers.  The  differences 
seen  between  the  Kipp  &  Zonen  and  the  Eppley  raise  doubts  about  their  absolute  accuracies  on 
time  scales  of  hours,  although  their  long-term  averages  are  quite  comparable. 
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1.  Introduction 


The  Upper  Ocean  Processes  (UOP)  group  at  the  Woods  Hole  Oceanographic  Institution  (WHOI) 
has  measured  a  suite  of  meteorological  parameters  from  buoys  and  ships  since  it  designed  the 
first  IMET  (Improved  Meteorology)  systems  (Hosom  et  al.,  1995).  The  suite  has  included  all  the 
parameters  required  for  making  estimates  of  air-sea  heat  and  momentum  fluxes.  Among  these 
parameters  is  the  downwelling  longwave  energy  flux.  Measurement  of  this  flux  is  challenging. 
Its  measurement  to  sufficient  accuracy  has  encountered  many  problems  and  generated  a 
considerable  literature  (Bums  et  al.,  2003;  Pascal  and  Josey,  2000;  Ji  and  Tsay,  2000;  Payne  and 
Anderson,  1999;  Fairall  et  al„  1998;  Philipona  et  al.,  1995;  Dickey  et  al.,  1994;  Miskolczi  and 
Guzzi,  1993;  Alados-Arboledas,  1988;  Albrecht  et  al.,  1974).  High  accuracy  measurements  at 
sea  have  increased  difficulties  over  land-based  measurements.  Recent  climate  research  programs 
have  required  increased  accuracy  of  air-sea  heat  flux  estimates.  The  UOP  group  has  striven  to 
steadily  improve  its  measurements,  and  the  work  described  here  is  part  of  this  attempt. 

Until  recently,  the  Eppley  Precision  Infrared  Radiometer  (PIR)  was  the  only  instrument  available 
for  measuring  downwelling  longwave  irradiance.  As  manufactured  by  Eppley,  the  PIR  is  not 
well  adapted  for  deployment  at  sea.  Its  aluminum  radiation  shield  is  fragile,  the  combination  of 
stainless  steel,  aluminum,  and  bronze  used  in  its  construction  do  not  stand  up  well  in  a  salt 
atmosphere,  and  it  is  difficult  to  mount  on  a  buoy.  Payne  and  Anderson  (1999)  described  the 
IMET  PER,  their  redesign  of  the  Eppley  PER,  and  its  calibration.  It  is  made  of  machined 
aluminum,  hard  anodized  and  sealed,  and  coated  with  a  tough,  durable  white  paint.  The  EMET 
PIR  is  manufactured  to  UOP  specifications  by  The  Eppley  Laboratory  and  has  been  successfully 
deployed  on  ships  and  buoys  for  approximately  fifteen  years. 

During  that  time  it  has  performed  reliably  but  satellite  transmissions  on  buoys  and  the  general 
high  level  of  electrical  noise  on  ships  has  caused  interference  with  the  very  small  output  voltages 
from  the  thermopile.  To  avoid  the  thermopile  voltages,  we  have  constructed  and  tested  a 
prototype  pyrgeometer,  henceforth  the  LWT,  in  which  all  measurements  are  with  thermistors; 
there  is  no  thermopile.  Since  thermistors  have  been  made  with  sufficient  stabilities  for  at  least  a 
decade  and  the  means  exist  to  calibrate  them  to  high  accuracy,  this  seemed  a  sensible  idea.  We 
also  present  some  results  of  what  is  probably  the  first  long  term  comparison  of  two  commercial 
longwave  flux  sensors:  the  Eppley  PER  and  the  Kipp  &  Zonen  CG4.  The  CG4  has  been  available 
since  early  2000. 

2.  PIR  Operation 

As  delivered,  the  Eppley  PER  contains  a  battery  circuit,  which  delivers  a  voltage  approximately 
proportional  to  the  incident  longwave  flux.  Fairall  et  al.  (1998)  has  shown  that  a  more  accurate 
estimate  of  the  incident  flux  can  be  made  by  removing  the  battery  circuit  and  recording  the 
thermistor  resistances  and  thermopile  voltage  directly  as  has  been  done  in  the  work  described 
herein. 

The  fundamental  radiometer  calibration  equation  for  a  pyrgeometer  with  a  thermopile,  and 
thermistors  in  the  case  and  dome,  (from  equation  10  in  Fairall  et  al.,  1998)  is: 
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(1) 


h =< +(t,  -n)6-(i-s, )/>>/«.  +b<j(t; -r‘) 

where  H  is  the  total  longwave  flux,  cris  the  Stefan-Boltzmann  constant,  Ts  is  the  absolute 
temperature  of  the  black  surface  of  the  thermopile  cap  (detector),  Tc  is  the  case  temperature,  £0 
is  the  emissivity  of  the  paint  on  the  cap,  p  is  the  infrared  reflection  coefficient  for  the  silicon 
dome,  k  is  the  thermal  conductivity  of  the  thermopile  mount,  B  is  a  second  calibration  constant, 
and  Td  is  the  absolute  temperature  of  the  dome.  Also, 


(T,-T')=ar  (2) 

where,  a  is  the  thermopile  sensitivity  constant,  694  K-V'1  for  the  PIR  (J.  Hickey  1993,  personal 
communication),  and  V  is  the  thermopile  voltage  in  volts.  Since  the  coefficient  of  (r,  -Tc)  is  a 
constant,  we  can  rewrite  the  equation  as 

H  =  or/+A(r,-rt)+s<r(r1,-7-;)  p) 

1  2  3 

These  three  terms  are  in  order  of  decreasing  magnitude  under  most  conditions.  Term  2  in 
equation  3  can  be  loosely  interpreted  as  the  thermal  conduction  from  the  detector  to  the  PIR  body 
and  term  3  as  the  radiation  exchange  between  the  detector  and  the  dome.  Hach  equation  requires 
only  three  temperature  measurements,  Ts,  Tc,  and  Td.  A  can  be  expressed  in  terms  of  s0,  the 
fundamental  radiometer  sensitivity  constant  by, 


For  the  Eppley  PIR  so  was  computed  from  se,  the  Eppley  calibration  constant,  by  equation  19 
from  Fairall  et  al.  (1998). 


_ £« _ 

{l  -  Accsea\[  +  B(l  -  /)]Tc3  } 


(5) 


where  f  =  0.5,  B  =  3.5,  and  Tc=  22  °C. 


Kipp  &  Zonen  give  calibration  constants  for  the  Albrecht  et  al.  (1974)  radiometer 
equation: 


H  =  &TcA+V/s 


(6) 


where 


s  1  =[k/T  +  4e0Tc3(l  +  s/T)]a 


(7) 
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A  photograph  of  the  Kipp  &  Zonen  CG4,  the  Eppley  PIR  and  one  of  the  thermistor  pyrgeometers 
is  shown  in  figure  1  and  cross-sectional  sketches  of  each  in  figures  2-4.  Both  the  Kipp  &  Zonen 
and  the  Eppley  have  a  thermistor  to  measure  case  temperature  and  a  thermopile  to  measure  the 
temperature  difference  between  case  and  the  blackened  detector  plate  under  the  dome.  The 
Eppley  has  a  thermistor  cemented  to  the  inside  of  the  dome  to  measure  dome  temperature.  The 
Kipp  &  Zonen  dome  is  quite  flattened  relative  to  a  hemisphere  and  is  thermally  tightly  coupled 
to  the  case.  Kipp  &  Zonen  believe  that  their  dome  configuration  contributes  a  negligible  amount 
to  the  radiation  sensed  by  the  detector  plate.  Thus,  the  last  term  in  equation  1  is  neglected  for  the 
CG4. 

3.  A  New  Pyrgeometer 

Three  fundamental  changes  were  made  from  the  Eppley  PIR.  All  measurements  are  made  with 
thermistors:  There  is  no  thermopile.  The  whole  sensor  is  made  of  machined  aluminum.  One 
advantage  is  that  the  thermal  conductivity  is  higher  than  either  bronze  or  stainless  steel,  which 
should  lead  to  smaller  temperature  gradients  within  the  sensor.  As  a  consequence  of  this  and  the 
white  paint  protecting  the  outer  surface,  we  are  able  to  dispense  with  the  radiation  shield.  These 
shields  are  fragile  and  do  not  perform  well  on  buoys.  The  blackened  disk  under  the  dome  is  2.2 
cm  in  diameter  filling  63%  of  the  area  of  the  dome  base.  The  equivalent  surface  in  the  Eppley 
PIR  is  1.1  cm  in  diameter  filling  17%  of  the  area  of  the  dome  base. 

The  external  shell  of  the  LWT  is  machined  aluminum,  hard  anodized  and  sealed,  and  coated  with 
Ameron  Siloxane  PSX  700,  a  rugged,  two-part,  bright  white  paint  containing  both  epoxy  and 
acrylic  polyurethane.  The  sides  are  at  a  45-degree  angle  to  the  base  to  better  reflect  shortwave 
radiation.  The  silicon  dome  was  purchased  from  the  Eppley  Laboratory. 

The  detector  plate  is  supported  by  three  1.6  mm  diameter  stainless  steel  rods,  which  are 
cemented  into  the  aluminum  internal  plug.  The  plug  is  also  hard  anodized  and  attached  to  the 
shell  by  four  screws.  Two  versions  were  assembled.  The  detector  plate  in  one  was  a  disk  of  1.3 
mm  thick  machined  of  aluminum;  in  the  other  it  was  a  0.1  mm  thick  aluminum  disk  cut  from  a 
soft  drink  can.  Both  disks  were  coated  on  their  upper  surfaces  with  Parson’s  Black  optical 
lacquer  purchased  from  the  Eppley  Laboratory. 

The  LWT  was  designed  to  accommodate  up  to  seven  thermistors.  Of  the  three  thermistors  potted 
in  the  base  plug,  one  was  2.5  mm  below  its  top  surface  under  the  middle  of  the  detector  plate, 
one  was  2.5  mm  below  the  bottom  of  the  detector  support  rods,  and  one  was  at  a  location  in  the 
plug  away  from  other  holes.  One  thermistor  was  cemented  to  the  underside  of  the  detector  plate 
with  conducting  epoxy.  Up  to  three  thermistors  can  be  cemented  to  the  inside  of  the  dome. 
LWT001  had  three  thermistors  approximately  120  degrees  apart  and  about  30  degrees  above  the 
dome  base.  LWT002  had  a  single  thermistor  about  30  degrees  above  the  dome  base.  For  test 
purposes,  multiple  thermistors  were  installed  to  monitor  temperature  gradients  within  the 
pyrgeometers. 

Two  types  of  thermistors  were  used,  both  glass-encapsulated  for  high  stability.  Thermometries 
type  BR55KB303K  has  a  diameter  of  only  1.3  mm  and  was  used  for  minimum  physical  size  and 
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thermal  mass.  Our  extensive  experience  with  Thermometries  epoxy-encapsulated  thermistors  is 
that  they  drift  less  than  2  mK/year.  Glass  encapsulation  is  well  known  to  improve  stability  of 
thermistors.  YSI  type  46032  has  a  diameter  of  2.2  mm  and  YSI  states  that  the  typical 
thermometric  drift  of  this  type  is  less  than  10  mK  in  100  months.  The  Thermometries  was  used 
in  the  dome  of  both  sensors  and  under  the  detector  plate  of  LWT002.  The  YSI  thermistor  was 
used  everywhere  else.  Both  types  have  a  nominal  resistance  of  30  kQ.  at  20°C. 

In  parallel  with  the  pyrgeometer  development,  we  built  two  IMET  modules  designed  to  record 
data  from  up  to  seven  thermistors.  Each  thermistor  has  a  0.1%  Vishay  resistor  in  series.  Five 
volts  are  applied  to  all  seven  resistor/thermistor  pairs  resulting  in  a  current  of  30-90  microamp 
depending  on  the  thermistor  temperature.  The  channels  are  individually  calibrated  using  a  set  of 
Vishay  0.01%  resistors,  equivalent  to  about  2  mK  for  both  thermistor  types  used  in  the  LWTs. 
Once  per  minute,  sixteen  scans  of  the  seven  thermistor  voltages  and  the  5  volt  exciting  voltages 
are  made  in  1.05  s,  and  the  averages  of  the  16  samples  are  recorded  on  a  PCMCIA  card  with 
time  and  date.  Resolution  of  an  individual  sample  is  equivalent  to  4  mK/count.  Power  dissipated 
in  the  thermistors  is  less  than  200  pW.  YSI  gives  a  dissipation  constant  of  10  mW/°C  in  stirred 
oil.  Assuming  that  the  potted  thermistors  have  at  least  as  good  heat  sinking  as  stirred  oil,  the 
short  time  that  the  current  is  applied  will  result  in  self-heating  of  less  than  0.02  K.  Similar  self 
heating  occurs  during  the  thermistor  calibration.  Altogether,  then,  we  expect  an  accuracy  of 
about  0.01  K  in  our  temperature  measurements. 

4.  Accuracy 


In  order  to  make  an  estimate  of  the  measurement  accuracy,  it  will  be  assumed  that  the  errors  are 
all  independent.  The  uncertainty  can  be  written  as: 


(»)*- 


—sr 

dr.  s 


V  ( 


dH 

dT 


ST : 


dr.*' 


| '-SB 
dB 


\ 2 


(8) 


Or 

{SH )2  =  ((Ts  -  Tc  )SA)2  +  {(4or(l  +  B)T3  +  a)ST5  }2  +  (aSTc  f  +  (4BoT]STd  f  +  [a{T4  -  Td4  )SbJ 

(9) 

From  our  thermistor  calibrations,  the  uncertainty  in  Ts  and  Tc  is  0.01  K.  That  in  Td  is  probably 
larger  from  self-heating  since  the  dome  is  a  poorer  heat  sink  and  from  gradients  (Philipona  et  al., 
1995)  so  we  estimate  that  as  0.1  K.  The  variation  in  A  and  B  through  a  number  of  repeat 
calibrations  of  Eppley  PIRs  is  about  5%.  For  LWT001,  (Ts-Tc)  is  about  0.4  K,  Ts  and  Td  vary 
between  275  K  in  winter  to  300  K  in  summer,  A  is  200,  B  is  3.7,  and  a(Ts4-Td4)  is  of  order  2.  H 
varies  from  300  to  400  W  m'2.  Combining  these  we  get, 

(5H)2  =  16.0  +  0.07  +  4.0  +  0.04  +  0.14  =  20.3  (W) 


4 


or  an  uncertainty  in  H  of  4.5  W  m'2,  about  1.5  %  which  agrees  with  the  estimate  by  Fairall  et  al., 
(1998)  for  the  Eppley  PIR.  The  dominant  term  in  this  calculation  is  from  the  uncertainty  in  A. 

5.  Calibration 

The  thermistors  were  calibrated  in  a  water  bath  to  an  accuracy  of  ±0.005  K  using  a  National 
Institute  of  Science  and  Technology  traceable  model  SBE35  thermometer  made  by  Sea-Bird 
Electronics.  The  data  were  fitted  to  the  Steinhart-Hart  (Steinhart  and  Hart,  1968)  equation: 

Y  =  a  +  b  log  (R)  +  c(log(R))3  (11) 

where  R  is  the  thermistor  resistance  in  ohms,  the  log  is  to  base  e,  and  a,b,c  are  the  calibration 
constants.  T  is  computed  in  kelvins. 

We  had  anticipated  calibrating  the  thermistor  pyrgeometers  by  the  method  of  Payne  and 
Anderson  (1999)  in  which  the  A  constant  is  computed  at  times  when  the  detector  plate 
temperature  is  equal  to  that  of  the  dome  and  LW3  in  equation  3  is  zero.  However,  the  slower 
response  times  of  the  new  sensor’s  detector  plate  meant  that  LW2  and  LW3  passed  through  zero 
at  approximately  the  same  time.  Instead,  we  arbitrarily  assumed  that  the  Kipp  &  Zonen  values 
were  the  most  reliable  and  did  a  least  squares  fit  of  equation  3  to  the  Kipp  &  Zonen  CG4  total 
longwave  fluxes  for  the  entire  data  set,  determining  values  of  A  and  B  for  LWT001  and 
LWT002.  Table  1  shows  the  calibration  constants  for  the  two  thermistor  pyrgeometers  and  the 
Eppley  PER.  LWT001  contains  the  1.3  mm  disk,  LWT002  the  0.1  mm  disk.  The  differences  in 
the  A  constant  are  striking. 


Table  1:  Calibration  Constants 


Sensor 

A 

B 

Eppley  26982 

366.0 

3.0 

LWT001 

200.4 

3.7 

LWT002 

101.2 

4.4 

6.  DataSet 

UOP  maintains  a  radiometer  calibration  facility  on  the  roof  of  the  Clark  Laboratory  on  the 
Quissett  Campus  of  the  Woods  Hole  Oceanographic  Institution.  This  facility,  shown  in  figure  5, 
has  a  totally  unobstructed  view  of  the  sky.  It  is  designed  to  allow  careful  leveling  of  up  to  8 
IMET  radiometer  modules,  as  well  as  mounting  up  to  5  standard  radiometers  to  be  recorded  on  a 
Campbell  Scientific  CR7  data  logger.  During  all  of  2002,  a  set  of  three  longwave  radiometers 
and  a  single  shortwave  radiometer  were  mounted  on  this  facility.  Thermistors  were  measured  by 
applying  a  precision  voltage  to  the  thermistor  and  a  Vishay  fixed  resistor  in  series  with  it  and  by 
digitizing  the  voltage  across  the  thermistor.  These  voltages  and  the  thermopile  voltage  were 
sampled  once  per  3  s.  One-minute  averages  were  recorded  by  the  CR7.  The  longwave  sensors 
included  a  recently  purchased  Kipp  &  Zonen  CG4,  an  Eppley  PER,  and  an  IMET  PIR  (Payne  and 
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Anderson,  1999).  The  shortwave  radiometer  was  an  Eppley  PSP  manufactured  by  Eppley  for  the 
UOP  group  to  a  design  similar  to  that  of  the  IMET  PIR,  i.e.,  white-painted  machined  aluminum 
with  no  shield.  The  equation  and  calibration  constants  supplied  by  Kipp  &  Zonen  were  used  for 
the  CG4.  Their  equation  is,  effectively,  equation  3  with  no  LW3  term,  and  with  Tc  substituted  for 
Ts  in  LW1.  The  Eppley  Laboratory  calibration  constant  of  the  PIR  was  used  to  compute  its  A 
constant  as  described  in  section  3  although  the  PIR  thermistors  were  calibrated  in  our  laboratory 
to  the  same  accuracy  as  the  LWT  thermistors. 

The  claim  by  Kipp  &  Zonen  that  their  dome  design  obviates  the  necessity  of  measuring  dome 
temperature  seems  to  be  valid:  A  careful  perusal  of  the  entire  data  set  showed  no  cases  of 
elevated  CG4  readings  over  the  PIR  caused  by  solar  radiation. 

Both  thermistor  pyrgeometer  modules  were  deployed  on  the  roof  for  nearly  a  year,  from  January 
27  through  December  31,  2002,  except  for  two  gaps:  18  April  to  30  May  and  7-29  August.  The 
thermistor  resistances  were  converted  to  temperatures  and  then  to  longwave  fluxes.  These 
module  files  and  the  processed  CR7  file  were  averaged  to  a  10-minute  series  for  analysis. 

7.  Anomalous  Data 

There  were  40  occasions  during  the  year  when  both  the  thermistor  pyrgeometers  recorded  total 
longwave  flux  values  that  were  20  W  m'2  or  more  higher  than  the  Kipp  &  Zonen  for  periods  of 
2-14  hours.  Most  of  these  periods  began  during  the  night  and  ended  sometime  after  sunrise. 
During  11  of  these  periods  the  IMET  PIR  also  showed  similarly  elevated  values,  although  of 
lesser  amplitude.  Since  most  of  the  anomalous  data  occurs  at  night,  we  ruled  out  shortwave 
leakage  of  the  domes  as  a  possible  cause.  Figure  6  shows  one  example  of  anomalous  LWT 
values.  In  Figure  7  both  the  LWTs  and  the  Eppley  PIR  are  anomalous  compared  to  the  Kipp  & 
Zonen  CG4.  There  was  no  rain  during  any  of  the  days  in  the  figures.  The  anomalies  in  each  case 
occur  in  LW2,  the  thermopile  voltage  term.  When  these  anomalies  occur,  there  is  no  identifiable 
common  feature  linking  those  instruments,  which  deviate  from  the  rest.  The  anomalous  data 
represent  6.9%  of  the  total  and  are  fairly  evenly  distributed  through  the  year. 

We  have  not  seen  another  long-term  comparison  between  the  CG4  and  the  PIR.  It  is  disturbing 
that  during  the  9.2  month  data  series  there  were  a  total  of  11  events  when  the  PER  values  were 
anomalously  high  by  20  W  m'2  or  more  compared  to  the  CG4.  In  each  case,  the  anomaly 
appeared  to  be  in  LW2,  the  term  containing  the  thermopile  voltage. 

8.  Comparison 


Since  the  Kipp  &  Zonen  does  not  have  a  term  analogous  to  term  3  in  equation  3,  and  since  it  was 
used  to  calibrate  the  LWTs,  we  compared  the  LWTs  with  the  Eppley  PIR.  Figures  8  and  9  are 
plots  of  LWT001  and  LWT002  total  longwave  flux  vs.  the  Eppley  PIR.  The  line  in  each  plot 
represents  a  least  squares  linear  fit  of  the  data.  The  agreement  is  within  approximately  5  W  m'2 
over  the  whole  range.  LWT002  has  somewhat  less  scatter  than  LWT001.  This  may  be  from  the 
thinner  detector  plate  and  smaller  thermistor  in  LWT002  yielding  a  slightly  shorter  time 
constant.  Table  2  gives  the  mean  and  standard  deviation  of  the  differences  between  the  two 
LWTs,  the  Eppley  and  the  Kipp  &  Zonen  for  the  anomaly-free  data  series  and  for  the  full  data 
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series  containing  the  anomalies.  Although  the  mean  differences  are  only  slightly  different,  the 
standard  deviations  are  nearly  doubled  when  the  anomalies  are  included. 


Table  2:  Pyrgeometer  mean 


differences  and  standard  deviations,  9.2  months  of  data 
Units  are  W  m"2 


No  Anomalies 

All  data 

Difference 

LWdiff 

SD 

LWdiff 

SD 

LWT001-K&Z 

-0.5 

4.4 

1.1 

8.5 

LWT001 -Eppley 

0.8 

5.3 

1.2 

8.2 

LWT002-K&Z 

0.9 

2.7 

2.3 

7.3 

LWT002-Eppley 

2.1 

3.1 

2.4 

6.7 

Eppley-K&Z 

-1.3 

2.5 

-0.1 

3.5 

It  is  encouraging  that  the  anomalies  do  not  affect  the  long-term  averages  appreciably  and  that  the 
long  term  means  agree  within  1-2  W  m‘2.  The  combination  of  these  mean  results  and  the 
anomalies  mentioned  previously  indicates  to  us,  however,  that,  although  long-term  averages  can 
be  accurate,  there  is  a  possibility  of  errors  amounting  to  tens  of  W  m'2  in  measurements  on  time 
scales  of  hours. 

We  also  carried  out  a  comparison  over  a  20-day  period,  which  contained  no  anomalies,  from  7  to 
26  October.  Figure  10  shows  the  total  longwave  flux  for  all  four  sensors,  the  CG4,  the  PIR,  and 
the  two  LWTs.  Figure  1 1  shows  the  difference  between  the  Kipp  &  Zonen  and  the  two  LWTs. 
Figure  12  shows  the  same  for  the  Kipp  &  Zonen  and  the  Eppley.  Figure  13  shows  the  shortwave 
flux  for  the  period.  Table  3  shows  the  mean  differences  and  standard  deviations  for  these  data. 
Differences  and  standard  deviations  are  of  the  same  order  as  for  the  whole  data  set. 

Table  3:  Pyrgeometer  mean  differences  and  standard  deviations,  year  day  280-299 

Units  are  W  m'2 


Difference 

LWdiff 

SD 

LWT001-K&Z 

-0.6 

4.2 

LWT001 -Eppley 

0.3 

4.4 

LWT002-K&Z 

1.1 

2.3 

LWT002-Eppley 

2.0 

1.5 

Eppley-K&Z 

-0.9 

1.7 

Lastly,  to  show  the  degree  of  agreement  between  the  four  pyrgeometers  and  the  effects  of 
ignoring  the  anomalies.  Table  4  contains  flux  values  averaged  over  the  whole  9.2  months,  with 
and  without  the  periods  containing  anomalous  data.  The  agreement  is  excellent  and  the 
anomalies  cause  an  elevation  of  the  PIR  and  LWT  averages  by  only  about  1  W  m'2  relative  to  the 
CG4. 
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Table  4:  Overall  longwave  flux  averages 
Units  are  W  m'2 


Sensor 

No  Anomalies 

All  Data 

Kipp  &  Zonen  CG4 

313.6 

312.8 

Eppley  PIR 

312.3 

312.7 

LWT001 

313.1 

314.0 

LWT002 

314.5 

315.1 

9.  Conclusion 

A  pyrgeometer,  intended  particularly  for  deployment  on  ships  and  buoys,  has  been  designed,  and 
two  have  been  built  and  tested  in  a  comparison  with  a  Kipp  &  Zonen  CG4  and  an  Eppley  PIR. 
The  comparison  data  set  consists  of  a  total  of  9.2  months  of  1-minute  averages  recorded  between 
January  and  December  2002  in  Woods  Hole,  Massachusetts. 

When  averaged  over  the  whole  data  set,  all  four  pyrgeometers  are  consistent  to  2.2  W  m'2  or 
better.  Standard  deviations  of  the  differences  of  the  two  LWTs  relative  to  the  Kipp  &  Zonen  are, 
however,  at  least  twice  that  of  the  Eppley.  The  standard  deviations  of  all  three  represent 
significant  differences  on  time  scales  of  hours.  The  LWTs  have  many  more  periods  of  such 
anomalies  than  does  the  Eppley.  Since  nearly  all  of  the  anomalous  periods  begin  before  dawn, 
shortwave  radiation  cannot  be  involved.  Several  of  the  periods  occurred  during  a  very  dry, 
rainless  period  so  moisture  on  the  outside  of  the  domes  is  unlikely  to  be  a  factor.  We  intend  to 
continue  our  observations  by  recording  additional  environmental  parameters  such  as  relative 
humidity,  precipitation,  and  wind  speed  to  try  to  find  the  cause  of  the  anomalies. 

We  conclude  that  contemporary  thermistors  allow  temperature  measurements  of  sufficient 
accuracy  that  the  thermopile  can  be  eliminated  from  pyrgeometers  although  our  design  is  not 
optimum.  We  also  conclude  that  differences  seen  between  the  Kipp  &  Zonen  CG4  and  the 
Eppley  PIR  raise  doubts  about  their  absolute  accuracies  on  time  scales  of  hours  although  their 
long-term  averages  are  quite  comparable. 
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Figure  2:  Cross-section  sketch  of  the  all-thermistor  pyrgeometer. 
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Figure  3:  Sketch  of  the  Kipp  &  Zonen  CG4  pyrgeometer.  Labeled  parts  are:  (A)  desiccant;  (B) 
white  sun  screen;  (C)  body  temperature  thermistor,  (D)  thermal  detector;  (E)  silicon  window;  (F) 
sub  frame. 
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Figure  4:  Sketch  of  the  Eppley  PIR.  (A)  top  and  dome  mount  (stainless  steel);  (B)  pyrgeometer 
body  (cast  bronze);  (C)  hole  in  which  body  thermistor  is  potted;  (D)  thermopile;  (E)  thermal 
detector;  (F)  internal  radiation  shield  (aluminum);  and  (G)  silicon  window. 
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Figure  5:  Radiometer  calibration  facility  on  roof  of  Clark  Laboratory,  Woods  Hole,  MA. 
The  Kipp  &  Zonen  CG4,  Epply  PIR,  Eppley/IMET  PSP  and  Eppley/IMET  PIR  are  in  the 
foreground.  An  LWT  is  at  the  top  of  each  of  the  two  long  tubes  just  beyond. 
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Figure  6:  Example  of  elevated  values  of  total  longwave  flux  of  LWT001,  LWT002  and 
IMET  PIR  (dotted  lines)  relative  to  the  Kipp  &  Zonen  CG4  (solid  line)  and  Eppley  PIR 
(dashed  line). 
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Figure  7:  Example  of  elevated  values  of  total  longwave  flux  of  the  LWT001 ,  LWT002 
and  the  IMET  PIR  (dotted  lines)  and  the  Eppley  PIR  (dashed  line)  relative  to  the  Kipp  & 
Zonen  CG4  (solid  line). 
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Figure  8:  Total  longwave  flux  of  thermistor  pyrgeometer  LWT001  vs.  Eppley  PIR. 
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Figure  10:  Total  longwave  flux  of  Kipp  &  Zonen  CG4  (solid  line),  Eppley  PIR  (dashed 
line),  and  LWT001,  LWT002  (dotted  lines),  7-26  October  2002. 
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Figure  1 1 :  Difference  between  LWT001  (solid  line),  LWT002  (dotted  line)  and  Kipp  & 
Zonen  CG4  total  longwave  fluxes,  7-26  October  2002. 
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Figure  12:  Difference  between  Eppley  PIR  and  Kipp  &  Zonen  CG4  total  longwave 
fluxes,  7-26  October  2002. 
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